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We investigate universality of the Nf = 2 running coupling in the Schodinger functional scheme, 
by calculating the step scaling function in lattice QCD with the renormalization group (RG) im- 
proved gauge action at both weak(M = 0.9796) and strong(M = 3.3340) couplings, where u = 
with ^SF being the running coupling in this scheme. In our main calculations, we use the tree- 
level value for 0{a) improvement coefficients of boundary gauge fields. In addition we employ 
the 1-loop value for them in order to see how scaling behaviours are affected by them. In the 
continuum limit, the step scaling function obtained from the RG improved gauge actions agrees 
with the previous result obtained from the plaquette action within errors at both couplings, though 
errors of our result are larger Combined fits using all data with the RG improved action as well 
as the plaquette action reduce errors in the continuum limit by 2% at the weak coupling and 22% 
at the strong coupling. 



The XXVI International Symposium on Lattice Field Theory 
July 14-19 2008 
Williamsburg, Virginia, USA 



* Speaker. 

^Riken BNL Research Center, Brookhaven National Laboratory, Upton, NY 11973, USA 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http://pos.sissa.it/ 



Universality of the Nf = 2 Running Coupling 



Keiko Murano 



1. Introduction 

The Nf = 2 p function in the SF scheme has been calculated recently with the plaquette action 
and the 0{a) improved fermion action in Ref.[|]] and it is found that the running coupling for 
Nf = 2 QCD becomes stronger than that for the pure gauge theory as the energy scale decreases. 
This behaviour is opposed to the perturbative prediction and thus is a genuine non-perturbative 
effect. In such a non-perturbative region, however, it is well-known that the calculation of the 
running coupling in the SF scheme becomes difficult, since the secondary minimum of the action 
comes close to the true minimum so that the auto-correlation time tends to be longer. Therefore it 
is important to check the non-perturbative behaviour of the j8 function mentioned above by using 
different gauge actions. 

In this study, employing the renormalization group (RG) improved gauge action, we have 
calculated the step scaling function in the SF scheme at both weak (m = 0.9793) and strong (u = 
3.3340) coupling regions, where u = with f sf being the running coupling in the SF scheme. 



2. Setup 



Using the similar setup given in Ref.|[I]], we consider the improved gauge action on an x T 
lattice in the SF scheme , 
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for a plaquette P^y and an 1 x 2 rectangular with the constraint that cq + 8ci = 1. For the 

RG improved gauge action, we take co = —0.331. To removed 0{a) scaling violations caused 
by boundaries, we have to tune the 0{a) improvement coefficients cf and cf , which have been 
calculated perturbatively at l-loop[Q] as 



CqC, 



cicfisl) 



c,{3/2 + cf^'^gl + 0{gt,)). 



(2.2) 



The scaling study of the SF running coupling for the SU(3) pure gauge theory with this gauge 
action [|3|], however, showed that the scaling violation at the strong coupling becomes larger for the 
1-loop value of cf ^ than that for the tree-level value. In this study we therefore take the tree-level 
value, cf = 1 and cf = 3/2. We have also performed an additional set of simulations with the 

PR PR 

1-loop value of c, ' , in order to see how scaling behaviour are affected by the choice of c, ' . For 
quarks, we employ the 0{a) improved Wilson fermion action (clover action) in the SF scheme[|I|], 
with the non-perturbative value of the improvement coefficient Qn^lQl and the 1-loop value of the 
improvement coefficient c,[Q]. 
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3. Simulation details 



We have calculated the step scaling function (SSF) in the weak coupling region (m = 0.9793) 
and the strong coupling region (u = 3.3340). Following the calculation procedure and the analysis 
in Ref.[|l|], we have calculated both g^p{L,a/L) and g\^{2L,a/L) at the same j8, in order to obtain 
the lattice SSF as 



Z{u,a/L) = glp{2L,a/L), u = gl^{L,a/L). 



(3.1) 



In order to make the continuum extrapolation as lim l.(u,a/L) = g(u), we repeat this procedure 

by changing L/a and j8 while keeping u = g|p(L,a/L) fixed. Throughout calculations, we tune the 
hopping parameter k: so that the PCAC quark mass m(xo) at xq = T /2 becomes zero at given L/a 
and ^0, where 



m{xo) 



I (dp + ^0 )/a (xq) + CAad^dofpixo) 
2fp{xo) 



(3.2) 



Throughout this study we take T = L. 

Errors of the SSF £ due to the small deviation of u from 0.9793 or 3.3340 are perturbatively 
corrected as 



'E,{u,a/L) = 'E,{u,a/L) + 1,' (u,a/L) x (m — u). 
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Similarly, errors of the SSF due to the small non-zero value of the PCAC quark mass are corrected 
as 



ll(u,a/L,z) = T.(u,a/L,0) + -:r-L(u,a/L,z) 
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(3.6) 



where the right-hand side of eq.(|3.6[) is taken from [|6|]. In eq( [3.5| , |3.6D , we make a z dependence of 
£ explicit, though hereafter we denote it as L{u,a/L) instead of L{u,a/L,z) for simplicity. 

We have carried out computations of the SSF at L/a = 4, 6, 8 and 2L/a = 8, 12, 16. Details of 
simulation parameters and analysis will be given in Ref.[^. Numerical simulations are performed 
on a cluster machine "kaede" at Academic Computing & Communications Center, University of 
Tsukuba. 



4. Some remarks for numerical simulations 

In the previous study with the plaquette action[|8], [l]], it has been reported that the auto- 
correlation time of the HMC updating tends to be longer in the strong coupling region. This long 
auto-correlation seems to be caused by rare but large fluctuations of the gauge part of 1 /g|p, which 
appear as a result of the tunneling between the true minimum and the secondary minimum. In 
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order to make the auto-correlation shorter, the modified gauge force which enhances a rate of such 
tunnelings has been introduced together with the reweighting method |l|]. 

We have checked whether a similar problem exists in the case of the RG improved action, by 
examining distributions of the gauge part of Fortunately we do not observe such rare but 

large fluctuations in distributions, and an example of distributions is shown in Fig|T[ We however 
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Figure 1: A distribution of the gauge part of l/g|p with the RG improved action, in the case of L/a = 16, 
j5 =2.755 and K- = 0.13334. 

observe that the distribution of 1 /g|p with the RG improved action has much wider width than that 
with the plaquette action[^. The wider width indicates that, while the HMC with the RG improved 
action samples configurations including ones near the secondary minimum more effectively, sta- 
tistical fluctuations of l/f|p with the RG action become also larger than those with the plaquette 
action. 



5. Result 

In Fig.§, we compare scaling behaviours between the (naive) lattice SSF L{u,a/L) and the 1- 
loop improved one 'L^^\u,a/L) for the RG-improved gauge action with the tree-level value of c, at 
the weak coupling (u = 0.9793, left) and the strong coupling (u = 3.3340, right). Here L^'^\u,a/L), 
whose lattice artifacts are perturbatively removed at k-loop, is defined by 

i:«f„a/n = L{2,u,a/L) 



where 5„ (n = 1,2 . . .) is given by 

'L{u,a/L) — (7{u,a/L) 



5i{a/L)u + &i{a/L)u^ + (5.2) 



o{u,a/L) 

and 8\{a/L) for the RG improved action is given in Ref.[^, |3|]. Since the scaling violation of 
is milder than that of Z, in particular at m = 0.9793, we hereafter use £*^') for the continuum 
extrapolation. 
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Figure 2: Scaling behaviours of E(open squares) and E'''(solid squares) at the weak coupUng (m — 0.9793, 
left) and at the strong coupling (u — 3.3340, right). 
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Figure 3: Linear continuum extrapolations of E'')(M,fl/L) at m = 0.9793 (left) and u ~ 3.3340 (right), 
together with the result of the ALPHA Collaboration[|l}, 



In Fig.^, we compare our Z^'^ with Z'^^ from the ALPHA Collaboration, where Z'^) is the 
2-Ioop improved lattice SSF, and observe that our errors of Z'^^ are larger in general. We think that 
these larger errors are caused partly by lower statistics and partly by wider widths of distributions 
mentioned in the previous section. In the figure we also plot linear continuum extrapolations of 
Z(i), which give a(M = 0.9793) = 1.061(19) with ;t;Vdof = 1.74 and a(3.3340) =4.91(58) with 
X^/dof = 0.585. These values are consistent with results from the ALPHA Collaboration, a{u = 
0.9793) = 1.072(4) and a(3.3340) = 5.60(16), though our errors are much lager as expected from 
errors of Z^^^. Note however that results from the ALPHA Collaboration are obtained from the 
combined fit with Z^^^ at all u by the form o{u) + p^^'°°PM'^(a/L)^ neglecting a tiny a/L term, 
therefore errors tend to be smaller than those form the individual fit. Here X is the order of the 
improvement coefficient c^: 1-loop for data at u in the weak coupling region and 2-loop in the 
strong coupling region[[T]], and the scaling violation starts at (3-loop) in lP-\u)/u. 

In order to examine scaling behaviours of Z^^' with the RG improved action more precisely, 
we have repeated the calculation of the SSF using the 1-loop value of boundary improvement 
coefficients cf'^ with condition B[|2|]. Z^'^ from the 1-loop cf'^ is plotted in Fig.Q together with 
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Figure 4: Combined fits to data including the tree c, (squire) at all L/a and the 1-loop q (cross) except 
L/a = 4, at u = 0.9793(left) and at u = 3.3340(right). Results from the ALPHA Collaboration ai-e also 
plotted(open circles). 



results from the tree c, and from the ALPHA Collaboration||i|]. We observe that three results show 
mild scaling violations for all a/Ls' at the weak coupling(the left figure) , while, in the strong 

PR 

coupling region (the right figure), the scale violation of the result from the 1-loop becomes 
too large at L/a = 4 to make the reliable continuum extrapolation with this point. A possible 
reason for this large scaling violation is that the perturbative estimate of c, ' becomes unreliable 
in the strong coupling region of the RG-improved action due to the large value of the bare gauge 
coupling ^Q. The p value which gives the same u is much smaller for the RG-improved action 
than for the plaquette action. For example, at m = 3.3340 and L/a = 4, the RG-improved action 
needs /3 = 2.1361 , which corresponds to gg = 2.8089 and gives a large(10%) 1-loop correction 
to cf, cf^^^gQ = 0.146. On the other hand, the tree-level value, cf =1.5, remains unchanged for all 
j3. Since is defined through boundary observables, a small difference in cf'^ may have a large 
effect on it. 

PR 

We now perform a combined fit using all our data with tree-level and 1-loop values of , 
except one with the 1-loop cf^ at L/a = 4 and u = 3.334. At the weak coupling( u = 0.9793), we 
use linear functions in L/a for data with the tree-level value of cf''^( 3 data) and the 1-loop values of 

PR 

Ct ' (3 data), while at the strong coupling(M = 3.3340), we use a linear fit for data with the tree-level 
value of c,( 3 data) and a constant fit for data with the 1-loop value of q excluding one at L/a = 4( 
1 data). From the simultaneous fits shown in Fig. |^ we obtain a{u = 0.9793) = 1.057(16) with 
xVdof = 1.22 and a{u = 3.334) = 5.57(22) with z^/<iof = 2.02. While both central value and 
error are almost unchanged at the weak coupling, an agreement with the result from the ALPHA 
Collaboration becomes better with smaller error at the strong coupling. 

We finally make a more complicated combined fit, using our data and fitting functions men- 
tioned above and employing the form a{u) + p^^'°°PM^(a/L)^ as fitting functions to data at all u 
from the ALPHA Collaboration [jl]]. We obtain a{u = 0.9793) = 1.0724(43) and a{u = 3.334) = 
5.559(125) with ;^^/dof = 2.39 in the continuum limit. Compared with result from ALPHA Col- 
laboration, 1.0736(44) for weak coupling and 5.60(16) for strong coupling, these errors are reduced 
by 2% and 22%, respectively. 
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6. Conclusion 

We have calculated the step scaling function in the SF scheme at weak (u = 0.9793) and 
strong (u = 3.3340) coupling regions with the renormalization group (RG) improved gauge action. 
Extrapolated values of the SSF from the RG improved gauge action agree with those from the 
plaquette action within errors at both couplings, though errors of the former are larger. 

From a combine fit using all data including ones with the plaquette action from the ALPHA 
Collaboration!!!]], we obtain o{u = 0.9793) = 1.0724(43) and o{u = 3.334) = 5.559(125) in the 
continuum limit. These errors are reduced by 2% and 22% from previous results in Ref.[[I]]. 

Finally we make two comments on calculations of the SSF with the RG improved gauge 
action. Firstly, it is better to use the tree-level value of 0{a) improvement coefficients cf'^ than 
the 1-loop value, in particular, in the strong coupling region. Secondly, the HMC algorithm with 
the RG improved action samples configurations including ones near the secondary minimum better 
than that with the plaquette action, though statistical fluctuations of l/f|p with the RG improved 
action become larger. 
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